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We investigate the dynamics of single and multiple scalar fields with exponential poten- 
tials, leading to power-law and assisted inflation, in loop quantum cosmology. Unlike in 
the classical theory, dynamical trajectories in loop quantum cosmology are generically non- 
singular, with a big bounce replacing classical big bang in the Planck regime. Post bounce, 
fvi ' after a phase of super-inflation, dynamical trajectories evolve towards the classical attractor 

7-H ■ in the inflationary scenarios. Depending on the initial conditions, bounce is shown to occur in 

kinetic as well as potential dominated regimes. We analyze the number of e-foldings resulting 
from the phase of super-inflation, and find the dependence of the maximum possible number 
of c-foldings on the equation of state at the bounce and on the steepness of the potential. 
We find that if the potential is not steep, this phase can lead to large number of e-foldings 
in power-law inflation. For the assisted inflation scenario, an increase in the number of fields 
I/"") . can yield a significant increase in the number of e-foldings during super-inflation. 

cr 

Jh ' I. INTRODUCTION 

OJQ; 

The inflationary paradigm provides one of the simplest ways to describe various aspects of the 
physics of the early universe in standard cosmology. However, despite its successes some of the 

0^ ' fundamental questions yet remain to be answered in this paradigm. These include the problem 

of initial singularity and the past incompleteness of inflationary spacetimes, solution of which are 

£T) [ expected to be provided by a quantum theory of gravity. Another issue is to understand the nature 

~r\ ' of inflaton itself. Though a simple <f> 2 inflaton potential fits extremely well with the observations, 

the origin of such a potential from high energy physics remains unclear. Theories invoking higher 

dimensions and supersymmetry yield multiple fields in the effective 4-dimensional description of 

• • \ the early universe. Further, models of supergravity naturally lead to exponential potentials [l|], 

. £h ' which do not yield exponential inflation as in the (/^potential. Such potentials lead to a power-law 

inflation provided the potential is not too steep [2J. The limitation of steep potentials can be 
alleviated using multiple fields, which can cooperate via Hubble damping to yield assisted inflation 
[3]. Interestingly, dynamical evolution in assisted inflation can also be understood as an effective 
single field power-law inflation. 

The dynamical properties of the positive exponential potentials leading to inflation in the 
Friedmann-Robertson- Walker (FRW) model have been widely studied [7H19J]. 1 In the forward 
evolution in the expanding branch, for sufficiently flat potentials, exponential potentials yield a 
scaling solution which is a late time attractor [7H9|. However, in the past evolution, dynamical 
trajectories in the power-law and assisted inflation scenarios lead to a kinetic energy dominated 
regime as the big bang singularity at vanishing scale factor is approached. Higher dimensional 
and supergravity models which yield these exponential potentials do not provide a generic remedy 
to overcome this fundamental problem. It has been hoped that incorporation of non-perturbative 
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Negative exponential potentials also lead to a rich physics, such as in Ekpyrotic scenario [4j. For an analysis of 
such a potential on the lines of the present work in loop quantum cosmology, we refer the reader to Refs. [g, |6j. 



quantum gravitational effects would lead to the resolution of the singularity problem. New physics 
resulting from such a non-perturbative quantization can also potentially affect the inflationary 
paradigm itself, and in principle may alleviate the observational constraints on power-law and 
assisted inflationary models. 

In recent years, quantization of homogeneous models in loop quantum cosmology (LQC) have 
shown that the cosmological singularities in various models can be successfully resolved [20|. LQC 
is a canonical quantization of cosmological spacetimes based on loop quantum gravity (LQG), a 
background independent and non-perturbative approach to quantize gravity. A key prediction of 
LQG is that the continuous differential geometry of the classical theory is replaced by a discrete 
quantum geometry at the Planck scale. The discreteness of quantum geometry has various novel 
implications. The Hamiltonian constraint, whose vanishing determines the physical solutions, is 
a quantum difference operator in LQC, and the resulting evolution equation is non-singular with 
uniform steps in the physical volume of the universe. The step size of the difference equation is 
determined by the minimum allowed eigenvalue of the area operator in LQG. At small spacetime 
curvature, the quantum difference equation can be approximated as a differential equation, and 
the solutions of the Hamiltonian constraint in LQC approach those of the Wheeler-DeWitt theory. 
In contrast, significant departures between LQC and Wheeler-DeWitt theory occur when space- 
time curvature approaches Planck regime. A state peaked on the classical trajectory when evolved 
backwards using quantum Hamiltonian constraint does not encounter a big bang singularity, but 
bounces when the energy density reaches a critical value p = /0 ma x ~ 0.41ppi (where ppi denotes 
Planck density) in the isotropic models 2114231] . 2 The existence of the bounce has been demon- 
strated to be a robust feature of the theory [25J . The singularity resolution at the quantum level 
occurs independent of the initial conditions or the choice of parameters for various matter models, 



including models with a non- vanishing spatial curvature 261-4281] . cosmological constant 29|, l3ll |. 



massive scalar field [231 ] and anisotropic models 3314351] . Using consistent histories approach, it has 



also been shown that the probability that a loop quantum universe encounters a singularity is zero 

For states which lead to a macroscopic universe at late times, discrete quantum evolution in LQC 
can be approximated b y a n effective continuum spacetime description using geometric formulation 



of quantum mechanics 39J-441I] . Effective dynamics has been widely used in LQC to understand the 
physical implications in various models and has provided important insights on the genericity of 
singularity resolution in LQC. In particular, analysis of isotropic and anisotropic models indicate 
that all strong curvature singularities in LQC are resolved 42H44J]. Using effective dynamics, effects 
of the loop quantum geometric corrections on eft 2 inflationary model have also been studied [5j, |45j- 

f. The phase of inflation turns out to be an attractor of dynamical trajectories after bounce 
and is preceded by a phase of super-inflation, where the rate of change of the Hubble rate is 



positive, for all matter obeying the weak energy condition 4a]. These properties have also been 
used to understand the likelihood of the occurrence of inflation for the (ft 2 potential 5ll453l]. 

In this work we analyze the dynamics of exponential potential leading to power-law and assisted 
inflation scenarios. We show that independent of the choice of initial conditions, effective dynamics 
is devoid of any singularities and power-law and assisted inflation fits naturally in LQC. In LQC, 
the scale factor bounces when energy density becomes equal to /0 ma x 5 leading to a phase of super- 
inflation which exits to the phase of classical inflation or decelerated expansion, depending on the 



2 The occurrence of bounce at p max is only a property of isotropic cases in LQC. In anisotropic models, bounces of 
the anisotropic scale factors are characterized by the behavior of energy density of the matter field and the shear 
scalar [241 ]. Bounces in anisotropic models in LQC can occur even when energy density and shear scalar do not 
attain a maximum. 
In contrast to the result in the Wheeler-DeWitt theory, where the answer turns out to be unity [37I.I38J). 



steepness of the potential. Analysis of the phase space trajectories in <j> — <p plane shows that in 
the regime p < Pmax/2, dynamical trajectories in LQC converge towards a curve which approaches 
the late time attr actor in the classical theory in the forward time evolution when p <C p max . In the 
backward time evolution, dynamical trajectories in the pre-bounce regime, approach the classical 
attractor at very early times in the contracting branch. In the classical theory, since contracting and 
expanding branches are disjoint, classical solutions either approach the forward time attractor (in 
the expanding branch) or are repelled by the forward time repeller (in the contracting branch). In 
LQC, quantum geometric effects lead to a non-singular evolution between the above attractor and 
repeller. In the classical theory, in the backward evolution of the expanding branch, dynamical 
trajectories approach those of the massless scalar near the singularity (in agreement with the 
Foster's theorem [54|). Hence, classical dynamics is governed by kinetic energy of the scalar field 
near the big bang/crunch with the equation of state w approaching unity. In contrast, in LQC, 
we find that the equation of state with the positive exponential potential can range from w = 1 
to w = —1 at the bounce, depending on the initial conditions. In LQC, though Foster's theorem 
is not applicable, we find that it is easier to find trajectories with a bounce where kinetic energy 
dominates over potential energy, i.e. with equation of state w ~ 1 in the backward evolution of 
the expanding branch. 

Since LQC leads to a phase of super-inflation preceding the inflationary phase, there are e- 
foldings which have a pure quantum gravitational origin. In the 4> 2 potential, unless at the bounce 
equation of state is w ~ —1, the number of e- foldings in the super- inflationary regime turn out to be 



very small [50j, |52j|. A pertinent question is what specific conditions on the exponential potentials 
will lead to a significant number of e-foldings during super-inflation and whether these potentials 
can further assist in the assisted inflation scenario. We find that for the exponential potential, 
the number of these e-foldings depends in a non-linear way on the equation of state at the bounce 
and the steepness of potential. As the steepness of the potential decreases, the maximum possible 
number of e-foldings during super-inflation occur for w —> — 1. However, for a non-zero value 
of the steepness parameter for the exponential potential, the maximum number of e-foldings in 
super- inflation do not occur for w = —1. Further, number of e-foldings can be significant even if w 
is not approximately —1 at the bounce. We analyze the issue of number of e-foldings during super- 
inflation for multiple field scenario and find that the number of e-foldings during super-inflation 
almost always turns out to be smaller than the equivalent single field scenario. Our investigations 
also reveal that by increasing the number of fields in assisted inflation scenario, the number of 
e-foldings during super-inflation phase can increase. 

In the following we begin with an overview of derivation of classical Friedmann and Raychaud- 
huri equations using the connection and triad variables of loop quantum gravity. We show the way 
these equations can be obtained in the canonical framework. These are then used to derive the 
modified Friedmann and Raychaudhuri equations using effective Hamiltonian in Sec. IIB. In Sec. 
Ill, we use the effective dynamical equations to investigate the case of power-law inflation. A de- 
tailed analysis of the behavior of dynamical trajectories and super-inflationary phase is presented, 
along with the dependence of number of e-foldings during super-inflation on steepness parameter 
of exponential potential and the equation of state at the bounce. This is followed by the analysis 
of multiple fields in Sec. IV leading to assisted inflation. Since system with multiple fields can be 
expressed as an equivalent single field case in assisted inflation, various results obtained in Sec. Ill 
immediately apply to the cases in Sec. IV. We conclude with a summary of results in Sec. V. 



II. CLASSICAL AND EFFECTIVE DYNAMICS FOR k = ISOTROPIC MODEL 

In this section, we review the way dynamical equations in the classical theory and the effective 
spacetime description of LQC can be derived in the Hamiltonian framework. The spacetime metric 
for the spatially flat model with lapse N chosen as unity is given by 

ds 2 = -At 2 + a 2 {t)(Ax{ + Ax\ + dx|) , (2.1) 

where a(t) is the scale factor of the universe. We will consider the spatial manifold to be M 3 with 
a physical metric q a b- Since the spatial manifold is non-compact, in order to define symplectic 
structure we introduce a fiducial cell C, which plays the role of an infra-red regulator to control 
divergences in the spatial integrals. 4 The physical volume of C is V = a 3 V , where V is the volume 
of the fiducial cell with respect to the fiducial metric q a b, given by q a b = ct~ 2 Qab- 

The gravitational sector of the classical phase space in LQG is spanned by the Ashtekar-Barbero 
SU(2) connection A l a and the conjugate triad Ef. In LQC, due to the underlying symmetries of 
the isotropic spacetime, these variables can be symmetry reduced as 

4 = c V-^ul E? = p^q V- 2 / 3 el , (2.2) 

where ef and °uj l a are the fiducial triad and co-triads associated with the cartesian coordinates. The 
triad p satisfies a kinematical relation, \p\ = V a 2 , where the modulus signs arises due to two 
possible orientations of the triad. In the following, without any loss of generality, we choose the 
fiducial volume of the cell C to be unity and fix the orientation of the triad to be positive. The 
triad obeys the following Poisson bracket relation with connection c: 

{c,p} = ^— ■ (2.3) 

Here 7 is the Barbero-Immirzi parameter, with value 7 ~ 0.2375 fixed by the black hole thermo- 
dynamics in LQG. In the following, we first discuss the classical Hamiltonian constraint and the 
derivation of Friedmann and Raychaudhuri equations in the canonical framework. This is followed 
by a similar derivation for the effective dynamics in LQC. 

A. Classical dynamics in Ashtekar variables 

Dynamical equations of the classical theory can be derived starting from the Hamiltonian con- 
straint expressed in terms of c, p variables and using Hamilton's equations. The classical Hamilto- 
nian constraint is given by, 

^cl = ~v4^P 1/2 ° 2 + Wmatt « (2.4) 

where % m att corresponds to the matter Hamiltonian, which in the case of minimally coupled scalar 
fields cpi, each with a potential Vi(4>i), is given by 



H 



matt 



E^att W = E [\ ^^ +P 3/2 Vmj ■ (2.5) 



4 The choice of the fiducial cell must not affect the physical predictions resulting from the quantization. For a 
discussion of these issues in LQC, and the way unphysical quantizations can be ruled out, we refer the reader to 
Ref. HHHl. 



The energy density pi and pressure Pi of a field fa, are defined as pi = % m att(i)/.P an d Pi 
— dl-L matt (i)/dp 3 < 2 , which turn out to be, 



-% + Wi), Pi = \%-Vi{4>i). (2.6) 



" 2p a ,vr *" ' 2p' 
Using the Hamilton's equations for matter variables we obtain, 

fa = {&,«*} = % , hi = {p^n c i} = -p^d^v^i) • (2-7) 

The time derivative of 4>i, then yields the Klein-Gordon equation for the field fa: 

fa + 3Hfa + d <Pi V i (fa)=0 (2.8) 

where H denotes the Hubble rate, H = a/ a. Using eqs. (|2.7p in expressions of pi and Pi and using 
(|2.8|) . it is straightforward to deduce that the total energy density p and pressure P satisfy the 
conservation law 

p + 3H(p + P) = (2.9) 

where, 



We now obtain the Friedmann and Raychaudhuri equations using the Hamilton's equations for 
the gravitational variables c and p. Hamilton's equation for the triad yields, 

P = {p,K d } = ~~Q~ c Ucl = —Pi ( 2 - u ) 

which using p = a? leads to the relationship between the connection c and the time derivative of 
the scale factor: c = ja. Similarly, Hamilton's equation for c leads to, 

c = {c,H d } = ^l^-Tia = ~lHa - AnG^aP (2.12) 

3 op 2 

where in the last step we have used the relation P = —dT-L ma ,tt/dp 3 ' 2 and the eq. (|2.1ip . Using the 
latter, the Friedmann equation follows directly by imposing the Hamiltonian constraint ()2.4p . and 
dividing it by the physical volume, 

c 2 8vrG<?W „ 2 8vrG 

7 z p 6 p a l A 6 

Similarly, using (|2.12p . we get the Raychaudhuri equation: 

a lc AmG „. , 

- = — = -— -p + 3P 2.14 

07a 6 

where in the last step we have used the expression of Hubble rate from fl2. 13f) . It is easy to 
verify that the Friedmann and Raychaudhuri equations together imply the conservation law ()2,9p . 
Thus, we obtain a closed set of dynamical equations for gravitational and matter variables. For 
the case of powerdaw and assisted inflationary scenarios, these equations break down as big bang 
the singularity at a — > is approached. In the regime close to the singularity, kinetic energy of 
the scalar field dominates the potential energy and solutions behave as in the case of a massless 
scalar field or equivalently the stiff fluid with equation of state w = P/p = 1. This behavior 
is in agreement with the Foster's theorem which shows that for positive potentials as considered 



here, the asymptotic behavior near the singularity mimics the case of a massless scalar 54]. As 
we will discuss below, the situation completely changes in LQC, where the big bang singularity is 
generically resolved due to underlying quantum geometric effects. 



B. Effective dynamics in LQC 

At the quantum level, evolution in LQC is determined by a quantum difference equation 
with a uniform step in volume. Using geometric formulation of quantum theory, where one 
treats the Hilbert space as a quantum phase space, it is possible to derive an effective Hamil- 
tonian for coherent states under controlled approximations 391441]. Extensive numerical simula- 



tions have shown that the dynamical trajectories resulting from effective Hamiltonian turn out 
to be in excellent agreement with the underlying quantum evolution for various matter models 
[22l . l23l . |26| . I27I . l29l . l3ll . l32j . |57| |. These effective equations incorporate underlying non-perturbative 



quantum geometric effect in a continuum spacetime description. In comparison to the classical 
Hamiltonian constraint, the key difference occurs in the gravitational part of the constraint where, 
instead of c 2 , the effective Hamiltonian constraint contains trigonometric functions of connection. 
For the spatially flat isotropic model, the effective Hamiltonian constraint is given by 

o 

H eS = #grav + flmatt = ~ Q n 2 - 2 |p| 1/2 sin 2 (/2 C ) + iJ matt « (2.15) 



where -ffmatt is considered to be arising from a Fock quantization of the matter Hamiltonian (|2.5[) , 
and ft = X/p 1 ' 2 . Here A 2 = 4(3^7)^, with lp\ as the Planck length, determined by the quan- 
tum geometry [33J. Since the quantum geometric effects do not modify the matter Hamiltonian, 
Hamilton's equations for matter yield the same dynamical equations as ()2.7p . Hence, for the above 
effective Hamiltonian constraint in LQC, the scalar field satisfies the Klein-Gordon equation (J2.8J) 
and the conservation law (|2.9|) . 

Unlike the dynamical equations for the matter variables, dynamical equations for the gravita- 
tional variables turn out to be qualitatively different. Hamilton's equations for the triad yield, 

8vrG7 d v 1/2 

p = {p, H cS } = -L -H cS = V -^ sin(2/2c) . (2.16) 

3 oc 7/x 

Using this equation, with the vanishing of the effective Hamiltonian constraint 

. 9 . . 8ttG\ 
sm 2( /2c ) = __ /0; (2.17) 

it is easy to derive the modified Friedmann equation: 

H 2 = f = 8KG p ( 1 __P_\ with pmax = _?_. (2 . 18) 



4p 2 3 ' V p m ax7 ' " mdX 8vrG 7 2 A 2 

Similarly, the time derivative of connection is given by 

c = {c, H cS } = ^2 jL HeS = _J_ sin (2 Ac ) c - -A^ sin 2 (/i C y/ 2 - 4vrG 7 aP . (2.19) 

It is then straightforward to obtain the modified Raychaudhuri equation by using the above equa- 
tion with the time derivative of fj2. 16() : 

" 47TG p (l - 4-^-) - 4.GP (l - 2-B-) . (2.20) 



a 6 \ Pmax/ V P 

Combining modified Friedmann and Raychaudhuri equations, one obtains another useful equation, 
measuring the time derivative of the Hubble rate 

H = -4:irG(p + P)(l-2-?—\. (2.21) 

\ Pmax / 



As expected from the dynamical equations of the matter variables, the modified Friedmann and 
Raychaudhuri equations together imply the conservation law (|2.9p . as in the classical theory. 

We now note some important features of the modified dynamics in LQC. Unlike the classical 
theory, where energy density can grow without any bound, p in LQC is bounded above by a 
maximum value p m3iX which is given by p max ~ 0-41ppi, where p-p\ is the Planck density. The 
Hubble rate attains a maximum value, |-ff| m ax = 1/(27 A) which is reached at p = p ma , x /2. Note 
that both /0 max and |-ff| max diverge as A — >• 0. This implies that the classical singularities are 
recovered in the limit where A, the parameter which captures the underlying quantum discreteness 
in LQC, vanishes. The boundedness of the Hubble rate has novel implications. First, it can be 
shown that it directly leads to the resolution of all strong singularities in spatially flat isotropic 
LQC |42| (see also Refs. 43|, [44| for the inclusion of spatial curvature and anisotropy). Secondly, 



there exists a phase of super-inflation, i.e. H > 0, when p ma , x < p < p max /2 14811 In literature, 
many interesting physical implications of this phase have been studied (see Ref. [20( for a review), 
including the qualitative dynamics of different potentials [5] and computation of probability for 
inflation in </> 2 potential to occur 52J, 153J]. 5 In the following, we will show the way this phase plays 



an important role in the qualitative dynamics for exponential potentials, which lead to power-law 
and assisted inflation scenarios. 



III. EXPONENTIAL POTENTIAL WITH A SCALAR FIELD 

In this section, we analyze the effective dynamics with an exponential potential. Such potentials 
arise naturally in various higher dimensional frameworks, such as in Kaluza-Klein theories and 
supergravity [1], and are of the form, 

V{4>) = v oe -^G^ (3-Q 

where V and k are constants. The exponential potential has various interesting properties in GR 
which have been widely studied in literature (see for eg. [2j, LZH10|) . The classical solution for the 
scale factor can be written as a power law [2j, [J, |8| : 

2 

a ex t n , with n = —^ . (3-2) 

Solving for the scalar field using classical equations, we obtain <j) oc (2/k) ln(t). The second time 
derivative of the scale factor implies that the power-law solution is inflationary in the classical 
theory if n > 1. It is non- inflationary if n < 1. Thus, in order to have an inflationary phase with 
the exponential potential, one requires the steepness parameter k to satisfy k 2 < 2. If the potential 
(|3.ip is steep, such that k 2 > 2, then it does not support inflation. Qualitative analysis with 
potential (13, ip shows that near the classical singularity, solutions are kinetic dominated and scale 
factor behaves as in the case of a massless scalar, i.e. a oc t 1 ' 3 [3, |l7|], as expected from the Foster's 
theorem [54j |. Further, for steepness parameter such that k 2 < 6, which constitute most interesting 
cases, there exists a potential-kinetic scaling solution which is a late time attractor [17|. 6 

In the classical theory, dynamical solutions for the potential (13.11) are generically singular. In 
contrast, in LQC, all the solutions with the above potential turn out to be non-singular. The 
differences between the classical theory and LQC only become significant in the Planck regime. In 



5 Dynamics of bounce and inflation has also been studied in models not incorporating LQC. See for eg. [581 ]. 

6 These properties also generalize to the case of the multi-field potential with n = YJ 2/kf (see for eg. [Ill.ll2l.ll8l]). 




FIG. 1: Evolution of the scale factor and energy density is shown in the Planck regime for LQC (solid 
curve) and the classical theory (dashed curve), for a single scalar field with potential (|3.ip . Values of the 
parameters of the potential are Vb = 0.001 and k = 0.5, and the units are chosen as Planck units. The scale 
factor in the classical theory vanishes, causing energy density to diverge. In LQC, energy density reaches a 
maximum value and the scale factor bounces. 
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FIG. 2: Long term evolution of the scale factor is shown for different initial conditions in the pre-bounce 
(left plot) and post-bounce (right plot). For the chosen initial conditions, bounce occurs at t — (not shown 
due to logarithmic scale). At very early (pre-bounce) and late (post-bounce) times, the classical power-law 
evolution is obtained. 



the regime p <C /9 m ax> the modified Friedmann and Raychaudhuri equations (12.1812. 20p approximate 
their classical counterparts (J2.13l2.14p . and, as expected, solutions in classical theory and LQC 
agree. This agreement occurs at late time after the bounce, as well as early times before the 
bounce, where LQC solutions furnish the same scaling property as in General Relativity (GR). A 
comparison of the solutions in the classical theory and LQC is illustrated in Fig. 1, and the early 
and late time behavior of solutions is shown in Fig. 2. In Fig. 1, we see that in the classical 
theory, the scale factor vanishes in a finite time in the past evolution of an expanding branch. 
This causes the energy density and the curvature invariants to diverge, and thus the physical 
evolution breaks down. In contrast, in LQC, starting with the initial conditions corresponding to 
the classical solutions at late times, we find that the dynamical solutions approximate the classical 
ones for a long time, however there are significant departures from the classical trajectory when 
energy density of the scalar field approaches p max . Using eqs. (l2.18|) and (I2.20p . we find that the 
scale factor in LQC bounces when p = p mSuX , and the backward evolution continues beyond the 
classical singularity. The linear relationship of In a versus lnt seen in Fig. 2 indicates the familiar 



power law behavior of the classical theory. Quantum gravitational effects encoded in eqs. (|2.18p 
and (|2.20p , thus lead to a past complete non-singular power-law inflationary model. These turn out 
to be robust features of the extensive numerical simulations performed with potential (|3.ip . In Fig. 
3, we show the plots for the relative errors in the effective Hamiltonian constraint (J2.15P for the 
solution discussed in Fig. 1. We see that the error in the effective Hamiltonian constraint is very 
small compared to the individual terms, H grSiV and -ffmatt- The behavior of relative errors shows 
that the effective Hamiltonian constraint is satisfied almost perfectly. (Relative errors remained 
similar in all the numerical solutions discussed in our analysis). 




FIG. 3: An analysis of the error for the solution with the solid curve in Fig. 1 is shown. The ratios 
|-£/eff/-Hgrav| and \H e g j H ma tt\ remain small during the evolution. 



Let us now examine the behavior of matter variables, <j> and (j), in the dynamical evolution. 
This can be understood by a phase space plot such as in Fig. 4, where we discuss the trajectories 
for the same values of parameters as in Fig. 1. Depending on the sign of the initial velocity 
(f>, classical curves approach the classical attractor in the expanding branch at large (ft at late 
times from the upper and the bottom halves of (j) — (j) plane. At early times, all of the classical 
dynamical trajectories, approach the regime where \(f)\ diverges, and near the singularity the phase 
space trajectories asymptotically behave as for the massless scalar model. Similar behavior of 
the trajectories exist for the contracting branch, where in the forward evolution classical curves 
starting from a repeller at very early times, approach the massless scalar regime near the big crunch 
singularity. In LQC, the behavior of trajectories in (p — (j) plane has some important differences with 
those in the classical theory. First, unlike in the classical theory, trajectories in LQC are bounded 
in |0| . This bound is a direct consequence of the upper bound on the allowed energy density p max in 
LQC. Second, in LQC, dynamical trajectories starting from the pre-bounce repeller in contracting 
branch approach the late time post-bounce attractor in the expanding branch. Using extensive 
numerical simulations, we found that dynamical trajectories in LQC first converge towards a non- 
classical curve in the post-bounce phase in the regime when energy density is p < /0 m ax/2, before 
they approach the classical attractor. This can be seen from the behavior of LQC trajectories in the 
upper half plane in Fig. 5, where small differences in the curve to which LQC trajectories initially 
converge and classical theory are visible when p ~ p m ax/2. The difference becomes negligible when 
the energy density decreases. A similar behavior exists for the repeller in the forward evolution of 
the contracting branch (in the lower half-plane of (f) — cp plot). 

We now consider some more detailed properties of the phase space trajectories in the Planck 
regime. In Fig. 5, initial conditions are provided, all corresponding to positive (j), in a narrow 
range of energy density in the super-inflation regime in LQC (p max > P > Pmax/2). We find that 
in the backward evolution, depending on the initial conditions on (ft and 4>, trajectories in LQC 
encounter a quantum bounce at p = p max , where cj> can be positive or negative. The magnitude of 
4> can take any value, including zero, at the bounce of the scale factor below the maximum bound 




FIG. 4: Phase plots with several initial conditions chosen in the 4> — <j) plane for a field with exponential 
function given by Vq = 0.001, k — 0.5. Initial conditions appear as solid dots on the solution curves, and 
arrows denote the direction of forward time evolution. The thick solid contour represents a level curve of 
energy density p = p max . The dashed contour represents energy density p = p max /2. 



4> 




FIG. 5: This plot shows a zoom of the phase space trajectories in the Planck region. The field and graphical 
notation is the same as in Fig. 3. Initial conditions correspond to densities between 0.6/9 max and 0.8/9 ma x- 
For the energy densities close to p m ax/2, LQC trajectories approach an attractor in the forward evolution, 
which approximates the classical attractor when energy density becomes much smaller than p m ax/2. 

set by p max . Thus, the quantum bounce can be kinetic (i.e. w = P/p ~ 1) or potential dominated 
(i.e. w ~ — 1). In the numerical simulations we found that if initial conditions are specified away 
from the bounce, trajectories which lead to kinetic dominated bounce are easier to find than those 
which lead to the potential dominated bounce. Thus in Fig. 5, trajectories approaching eft ~ at 
the bounce, in the backward evolution from the initial conditions specified in Planck regime, are 
fewer than those which correspond to (f) oo 0. We emphasize that this behavior strongly depends 
on the specification of the initial conditions in the phase space plot. Finally, we note that since in 
LQC (p an d </> can take on a range of values at the bounce, the equation of state for the scalar field 
at the bounce can take on values between w = 1 and w = — 1. 

Having noted the main features of the dynamics and the phase space trajectories, we now turn 
our attention to understand the super-inflation phase in this model. This phase occurs when energy 
density of the scalar field is /O max > p > /Omax/2, where H > 0, and contributes to additional e- 
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FIG. 6: The amount of super- inflation (in e- foldings) vs w, the equation of state, at the time of bounce 
for a field with steepness parameter k = 1. A wide range of w values is obtained by choosing many initial 
conditions along the curve defined by p — 0.95p max - 

foldings originating purely due to quantum gravity effects. A pertinent question is whether these 
e-foldings are large and under what conditions. It turns out that the number of e-foldings depend 
on the equation of state at the bounce, as well as the parameter fc in a non-trivial way. 7 For each 
given value of w with a particular value of k, there exist two values of the number of e-foldings in 
the super-inflationary regime, as is depicted in Fig. 6. In the special cases of w = —1 at the bounce, 
these two values coincide. This behavior can be understood by noting that for the trajectories in 
the (j) — (p plane, for a fixed value of (j> at the bounce, energy density can be equal to p max for two 
values of (j) with same magnitude but different signs. For different signs of (j) at the bounce, the 
evolution of energy density and Hubble rates in the regime p max > p > p ma , x /2 is different. This 
difference in evolution causes, two values of the number of e-foldings for the same equation of state, 
one corresponding to positive </> at the bounce, and the other to the negative (j). For a given value 
of w at the bounce, the magnitude of the difference depends on the value of the parameter k. 

This behavior is illustrated by the Figs. 7 & 8, where for a fixed value of the parameters in 
the potential, we have compared the phase space trajectory, the behavior of energy density and 
Hubble rate, and the amount of super-inflation versus the equation of state at bounce. This is 
done for two sets of initial conditions leading to the bounce at the same value of (j), and <p with 
equal magnitude but opposite signs. We find that in the case for <ft > at the bounce (Fig. 7), the 
number of e-foldings in the super-inflation regime is less than 0.5, whereas for the case when <p < 
(Fig. 8), it is approximately equal to 4. From the Fig. 7 & 8, we see that this difference is due to 
the relatively longer duration in which p > p max /2 when <ft is negative at the bounce, compared to 
when it is positive. The cause of the difference in this duration is ultimately tied to the fact that, 
in the forward evolution, trajectories for which the bounce occurs in the lower half of the (ft — <fi 
plane cover a longer path in the phase plane to reach p mauX /2 and merge with the forward attractor 
(which lies in the upper half plane) than those which bounce in the upper half plane. In Fig. 9, we 
illustrate the special case of w = — 1 at the bounce. In this case, the number of e-foldings during 
super-inflation turn out to be more than those for the initial conditions in Figs. 7 & 8, where the 
equation of state at the bounce is approximately zero. A comparison of the plots of p and H in 
Figs. 7, 8 &; 9, reveals that for the case of bounce with w = — 1, the variation of energy density 
and rate of change of the Hubble rate favors a longer duration of super-inflation and hence a larger 
number of e-foldings. 



These do not depend on the value of Vo, a change in which only causes a trivial shift in 
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FIG. 7: Phase space trajectory and the behavior of energy density and time derivative of Hubble rate, along 
with the number of e-foldings versus equation of state are shown for k = 0.16 and V = 0.1. A circle on the 
curve denotes the value at the bounce. Each dot in the plot of the number of super-inflationary c-foldings 
corresponds to an initial condition taken at p = 0.999p max . 




FIG. 8: For the same values of V and k as in Fig. 7, the phase space trajectory, behavior of energy density 
and H, and the value of equation of state at the bounce are shown for different initial conditions. Since 
4> < at the bounce, the number of e-foldings during super-inflation are higher than in Fig. 7. 

Interestingly, as the value of the parameter k is changed, the maximum number of e-foldings 
in the super-inflationary phase changes. Further, the value of the equation of state at the 
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FIG. 9: The plot shows the way variation in energy density and H, due to the behavior of 4> affects the 
number of e-foldings during super-inflation. The parameters of potential are same as in Fig. 7, and only 
initial conditions are changed. The bounce occurs at w = — 1, which yields more (but not the maximum) 
c- foldings in comparison to initial conditions in Fig. 7 & 8. 

bounce which yields the maximum number of super-inflationary e-foldings also depends on k, 
though as expected it is independent of V . From Fig. 6, we find that for the parameter k 
equal to unity, there are approximately 0.8 e-foldings as the maximum number of e-foldings 
due to super-inflation when the equation of state at the bounce is w ~ 0.5. Whereas, for 
k = 0.16 the maximum number of e-foldings are approximately 8 when w ~ —0.5 at bounce. We 
performed a detailed analysis of this behavior for various sets of initial conditions with different 
values of k, and the results are represented in Fig. 10. Two main features appear from this analysis: 

(i) The dependence of the maximum number of e-foldings in super-inflationary phase is non- 
linearly related to the value of k. As k decreases, i.e. the potential (|3.1j) becomes shallower, the 
maximum number of e-foldings during super-inflationary phase increases. In our simulations, we 
found that for k = 0.06, the number of e-foldings during super-inflation can be as high as 50. 
As the potential becomes steeper, i.e. k increases, the maximum number of e-foldings during 
super- inflation decreases. 

(ii) The equation of state at the bounce yielding the maximum number of e-foldings during super- 
inflation depends on A; in a non-linear way. For steep potentials, maximum number of e-foldings 
in super-inflation occur when the bounce is not potential dominated. As this parameter decreases, 
the potential becomes shallower, and the maximum number of e-foldings during super-inflation 
occur for equation of the state at the bounce approaching —1. Fig. 10 suggests that the maximum 
number of e-foldings occur for w = — 1 as k approaches zero. 
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FIG. 10: Left plot shows the dependence of the maximum number of e-folding during super-inflation (plotted 
along y-axis). This number increases as k decreases. The right plot shows the relationship between the 
equation of state at the bounce which produces the maximum number of e-foldings during super-inflation 
and the parameter k. 



IV. ASSISTED INFLATION 



As discussed earlier, for the potential of the form (|3.ip . if k 2 > 2, a single scalar field does 
not lead to inflationary dynamics in GR. However, multiple fields (pi, each with an exponential 
potential (|3.ip with a steepness parameter ki, which couple through Hubble expansion, can lead to 
inflation even if k 2 > 2 for each field. This novel mechanism, in which individual fields cooperate 
to produce inflation when none of them could have done so individually, forms the basis of assisted 
inflation [3j. The combined potential of multiple fields, V, is given by 



V = Y,Vi(<f>i) = Y,Vo< 



?S,iTGki(j>i 



(4.1) 



where Vj denotes the potential of each of the field (pi > and V is considered to be the same for all 
fields. In such a scenario, the dynamical equations for the gravitationally coupled multiple fields 
can be re-expressed as those of an effective single field (p with an equivalent k eq given by [3J] 



1 

i.2 



Z^ fr? • 



(4.2) 



Thus, even though each kf > 2, k 2 can be less than 2. If so, the multiple field configuration 
can successfully lead to inflation with a power law expansion a oc £™ cq with n eq = > 2/ kf > 1. 

i 

In this way different fields assist each other in causing inflation [3]. In fact, as the number of 
fields increases, inflation becomes possible for larger values of the individual fcj's, and the effective 
potential V becomes shallower. Since the multi-field system in assisted inflation reduces to an 
effective single field scenario, various results of the previous section are applicable here. In the 
classical dynamics, in the forward evolution the fields evolve with energy densities in fixed ratios, 
leading to a scaling solution which is an attractor [12|. In the backward evolution, as the past 
singularity is approached, dynamical trajectories in (p — <j> plane approach those of the massless 
scalar fields and the potential V becomes insignificant near the singularity. On the other hand, 
in LQC, singularity is avoided irrespective of the choice of the number of fields, parameters V 
and fcj's and the initial conditions. However, the behavior of each field near the bounce is now 
determined collectively by the steepness parameters ki and the initial conditions of all the fields. 

To understand the dynamics in a multi-field scenario, we performed various numerical simula- 
tions with a system of three fields. An example of the dynamical trajectories in the (p — (p plane 
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FIG. 11: The multi-field dynamics of 3 fields </>i, 02 and 03 is shown (ordered left to right). The fields have 
potential parameters as: V = 0.1, fci = 1.5, k% = 1.75, and k% = 2 respectively. The initial conditions are 
0i (0) = 0.4, 2 (0) = 0.7, 03 (0) = 1 and 0i(O) = 2 (O) = 3 (O) = 0.091 in Planck units, and their location is 
denoted by thick black dots on the curves. The arrows indicate the forward time direction, and the dashed 
lines represent classical solutions to the same initial conditions. The bounce occurs at the three points 
marked by an open circle. 




1.2 



FIG. 12: Plot of the dynamical trajectories in the classical theory (dashed curves) and LQC (solid curves) 
are shown for the three field system with choice of parameters as in Fig. 11. The initial conditions for 0i, 02 
and 03 are similar to those in Fig. 11, with the only different values being 02(0) = 0.08 and 03(0) = 0.7. (in 
Planck units). The trajectories for fields are plotted from left to right in the order: 02, 0i, 03. The thick 
dashed curve shows the dynamical trajectory in LQC for the field 03 treated as in single field evolution. 

are shown in Fig. 11 where we consider steepness parameters as k\ = 1.5, £;2 = 1.75 and k^ = 2.0. 
The equivalent steepness parameter for the single field case turns out to be k eq ~ 0.99. Thus, even 
though each of the three scalar fields can not individually lead to inflation, the effective steepness 
parameter is such that assisted inflation is possible. From Fig. 11 we see that in the backward evo- 
lution of the expanding branch, dynamical trajectories in LQC do not follow classical trajectories in 
to the singularity but bounce and turn-around to approach the attractor in the backward evolution 
of the contracting branch. In the forward evolution, dynamical trajectories in LQC approximate 
the classical trajectories when p -C /0 max and approach the late-time attractor of the multi-field 
system in the classical theory [121 ]. Unlike in LQC, in the classical theory, the attractor in the 
forward evolution of the expanding branch is disjoint from the repeller in the forward evolution of 
the contracting branch due to the presence of the big bang singularity. We find that the field 0i 
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approaches <f> = in the forward evolution more slowly than fa and fa. In terms of the single field 
evolution, this behavior can be understood by noting that the the effective steepness parameter 
k cq of the system is closer to k\ than &2 or £3 . 

The effect of the choice of initial conditions in multi-field scenario can be seen by comparing 
Fig. 12 with Fig. 11. In Fig. 12, we consider the same set of parameters and initial conditions 
as in Fig. 11 but with a change in the initial potential energy of the fields fa and fa. We find 
that the new initial conditions cause a non-trivial change in the way different fields approach the 
scaling solution at late times (after the bounce), and also at early times in the pre-bounce phase. 
Unlike the case in Fig. 11, the bounce of the scale factor occurs when fa is negative, and the 
kinetic energy of the field is less dominant than before, though the implications these changes have 
on the system is much more complicated than in the single field cases explored earlier. Fig. 12 also 
contrasts the behavior of the single field case with that of multi-fields. For the field fa, we show 
the dynamical trajectory (thick dashed curve) for the same initial conditions and parameters but 
evolved as an individual field. It is clearly seen that the attractor and the location of the bounce 
differ between the single field and multi-field cases with the same initial conditions, due to the 
additional contribution to Hubble damping by fields fa and fa in the latter case. 
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FIG. 13: The number of c-foldings obtained during super-inflation versus equation of state for a three field 
system in assisted inflation is shown. Parameters are k\ = 1.5, &2 = 1-6, k$ = 1.7 and V = 0.1. Points are 
shown for 2000 randomized initial conditions with randomized initial field densities. 

The next step in the analysis is to understand the effect of super-inflationary phase on the 
multiple field system. In particular, an important question is whether the number of e-foldings 
during super-inflation can further assist the assisted inflation scenario. We investigated the number 
of e-foldings allowed during super-inflation for potentials with various parameters in the case of a 3 
fields. An example of the results for steep potentials Vi is shown in Fig. 13, for k\ = 1.5, hi = 1.6, 
&3 = 1.7 and V = 0.1, where 2000 data points for number of e-foldings obtained by simulations 
with varying initial conditions are plotted. These were obtained as follows. First a pseudo-random 
initial density, p±, is chosen between 0.001p max and 0.99/5 max . Then p2 is chosen in a pseudo- 
random way with value between and p max — p\. Similarly, the energy density pj, is chosen. The 
densities are then permuted pseudo-randomly, to prevent p\ from being larger on average than 
the others. These are then paired with the fields <p,fa,fa, and a large array of initial conditions 
is generated along the density level curve defined by p = pi (where i = 1,2,3) on the portion of 
this curve nearest to the potential regime. This placement, along with the fact that the initial 
density generation method statistically favors high initial ptotai = Yl Pu an d helps produce bounces 
with total equation of state w 96 1. Each field is separately assigned an initial condition from the 
appropriate array generated in the previous step, and the LQC equations are solved with these 
initial conditions. These solutions are used to find the number of e-foldings during super-inflation, 

16 



and the process is repeated. The effective steepness parameter for the above case is k eq m 0.92. 
Comparing our result with the single field case (see Fig. 10), we find that for all initial conditions, 
the number of e-foldings during super-inflation is less than the maximum number of e-foldings for 
the single field case with the steepness parameter k = k eq . 
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FIG. 14: The left plot shows the number of e-foldings obtained during super-inflation for three fields with 
ki = 1.5. The right plot corresponds to 20 fields with the same fej's. In each plot 1000 initial conditions are 
considered, randomized as in Fig. 13 with the added constraint that the total initial density of each field is 
0.9p ma x, divided by the number of fields. 



For the three field scenario, though the number of e-foldings during super-inflation are not 
significant, they increase on addition of more fields. This is confirmed from the plots in Fig. 
14 where we have contrasted the three field system with 20 fields, each having same steepness 
parameter ki = 1.5. The initial density in both plots corresponds to 0.9/9 max , evenly distributed 
amongst fields and initial conditions randomized as above. Notable differences arise on changing 
the number of fields. First, the plots show that the amount of e-foldings during super-inflation 
can increase with an increase in the number of fields in the assisted inflation scenario in LQC. 
Second, we find that in the case of three fields, the number of e-foldings during super-inflation are 
closer to the value in the equivalent single field case (with k eq = 0.866), than for the 20 fields (with 
k eq = 0.335). Further, the distribution of points in the plot of 20 fields is in narrow region than in 
the case of three fields. Our simulations showed that this is due to lower initial energy densities of 
the individual fields for the case of 20 fields. 

Finally, let us understand the way the number of e-foldings during super-inflation change if one 
chooses a shallower potential. To examine this we performed various simulations involving ki < v2 
for three fields. An example of one of the results is plotted in the left figure of Fig. 14. Here the 
steepness parameters are equal for all the fields, k\ = &2 = ^3 = 0.3984 with k cq = 0.23. Following 
the process of providing initial conditions as for Fig. 13, e-foldings in the super-inflationary phase 
were obtained for 2500 numerical runs. We found that almost all points lie in the area defined by 
two curves for the single field cases, one with steepness parameter k = k eq = 0.23 and another with 
k' = ki = 0.3984. Further, data points display a tendency to clump near the boundaries of the 
limiting k and k' curves in some regions. This plot is to be contrasted with the second plot in Fig. 
14, obtained with different values of k\, &2 and k^ yielding k cq = 0.23. In this case the data points 
are obtained as follows. Note that the power in the asymptotic scale factor solution can be written 
as n eq = Y2 n i- This power can be randomized in the same way as the randomization of densities 
in Fig. 13. First, n\ is chosen pseudo-randomly between 0.001n eq and 0.99n eq . Then ni is chosen 
pseudo-randomly between and n eq — n\. Similarly, n% is chosen such that Yl n i = n eq- We then 
use the relation ki = y/2/m to determine the steepness of each field. Then the initial conditions 
are randomized as in Figure 13. This process was repeated 2500 times. We see that data points 
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are now found outside their previous boundaries in left plot in Fig. 14, on the bottom right and 
left, and are uniformly distributed in comparison. This behavior shows that the distribution of 
points in the plot of number of e-foldings versus equation of state is sensitive to the individual 
steepness of different potentials in assisted inflation. However, note that no points were found 
which led to more e-foldings than the maximum number of e-foldings during super-inflation for the 
single field case. One can thus conclude that maximum number of e-foldings in the single field case 
achieved during super-inflation provides a good upper bound on the number of e-foldings during 
super-inflation for the multi-field scenarios with exponential potential. 




FIG. 15: Left plot shows the number of e-foldings during super- inflation for three field system with k\ = 
&2 = &3 = 0.3984 and V = 0.1. The thick dotted upper curve corresponds to the e-foldings during super- 
inflation for the single field case for k = k cq — 0.23 and the lower one for k' = ki = 0.3984. The right plot 
shows the number of e-foldings during super-inflation for the three field system for V = 0.1 and fc C q = 0.23 
with differing k[s and initial conditions chosen in a pseudo-random way as described in the text. 



SUMMARY 



Scalar fields with exponential potentials appear naturally in various theories of gravity at high 
energies. These potentials have been used to construct models of power-law and assisted inflation. 
However, as in the case of the chaotic 4> 2 inflation model, the past evolution in these models is 
generically singular in GR. It has been hoped that incorporation of quantum gravity effects would 
resolve the initial singularity in these models and provide insights on the new physics. In this work, 
we probe this issue within the framework of LQC, which has been successful in resolving various 
cosmological singularities [20] . Using the effective spacetime description, we analyze the dynamical 
evolution of single and multiple fields with positive exponential self-interaction potentials. We 
show that generically non-singular power-law and assisted inflation scenarios can be constructed 
in LQC. Unlike the classical theory, where all solutions are singular, past evolution using modified 
dynamical equations in LQC generically result in a bounce of the scale factor when energy density 
becomes equal to p max ~ 0Alpp\. As an important element of new physics at the Planck scale, a 
phase of super-inflation which occurs when /? max > P > Pmax/2 precedes the classical inflationary 
phase. This phase contributes additional e-foldings of inflation to those in the classical inflationary 
phase, which can potentially affect the viability of these models by modifying the constraints on 
the steepness parameter in power-law and assisted inflation scenario. The number of e-foldings 
contributed by super-inflation can be large for non-kinetic dominated bounces. Our results show 
that as the steepness parameter k (or A; eq in assisted inflation) decreases, the number of possible e- 
foldings which can occur during super-inflation increases. In assisted inflation, this can be achieved 
by adding more fields without changing the steepness of the field potentials. This result indicates 
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that if the exponential potential is suitably modified to end the inflation, a given number of e- 
foldings can be achieved in LQC in principle with a lesser number of fields than in GR. Due to 
above reasons, quantum gravity effects as understood in LQC, can potentially further assist assisted 
inflation scenario. 

In the power-law and assisted inflation scenarios, at late times in the expanding branch, classical 
solutions approach an attractor, and at early times, in the backward evolution of the expanding 
branch, dynamical trajectories asymptotically approach those of the massless scalar field near the 
initial singularity. Dynamical trajectories in the contracting branch, disjoint from the expanding 
branch in the classical theory, in the backward evolution approach the attractor at early times. 
We find that in LQC, due to generic resolution of the big bang singularity, the situation changes 
dramatically. In the forward evolution of the expanding branch, dynamical trajectories in LQC 
approach the attractor of the classical theory. However, we find that before that occurs, dynamical 
trajectories approach a non-classical curve in the phase space for p ~ /9 m ax/2, which approaches 
the classical attractor as the energy density decreases. Existence of this curve is also found in 
the pre-bounce regime, in the backward evolution of the contracting branch, which approaches 
the backward time attractor at early times. Thus, quantum geometry effects in LQC bridge the 
pre-bounce repeller with the post-bounce attractor in the forward time evolution. 

In the classical theory dynamical trajectories asymptotically approach those of the massless 
scalar field near the initial singularity, and the approach to singularity is thus dominated by kinetic 
part of the matter Hamiltonian. In LQC, since all solutions bounce at a non-zero volume, solutions 
at the bounce can be dominated by kinetic as well as potential terms. Detailed analysis of the phase 
space trajectories reveals that the equation of state at the bounce varies from w = 1 to w = — 1 
depending on the initial conditions and the steepness of the exponential potential. The maximum 
number of e-foldings during super-inflation increases as the steepness parameter decreases. The 
dependence is non-linear, and the number of e-foldings can be very large for shallow potentials. 
An interesting feature of this model is that the maximum number of e-foldings during super- 
inflation, depending on steepness parameter, occur for different values of equation of state at 
the bounce. As the steepness parameter decreases, the equation of state at the bounce which 
yields maximum number of e- folding during super-inflation shifts towards w ~ — 1. For steeper 
potentials, the equation of state at the bounce yielding maximum number of e-foldings shifts 
towards w ~ 1. Analysis of the multi-field system in assisted inflation scenario shows that the 
number of e-foldings obtained during super-inflation are roughly bounded above by the single field 
case with an equivalent steepness parameter. We find that the likelihood of obtaining high e- 
foldings during super-inflation for a given set of initial conditions with same equivalent steepness 
parameter changes if the individual steepness parameters of multiple fields change. The number of 
e-foldings during super-inflation can increase with an increase in the number of fields in assisted 
inflation. 

An important question is the way above results are affected if we go beyond the homogeneity 
and isotropy assumption in loop quantum cosmology. The first step in this direction would be to 
analyze the dynamics with exponential potentials in anisotropic situations, such as in the Bianchi-I 
model. Due to availability of the underlying loop quantization 13311 . which also indicates a generic 
resolution of strong singularities at the effective spacetime level [441 ]. it is expected that the classical 
singularity will be avoided in a qualitatively similar manner as in the isotropic model. However, 
due to non-zero anisotropic shear, the resulting physics at the Planck scale is expected to be 
richer. As one of the features, the approach to the classical singularity would be affected by the 
presence of anisotropics. Unlike the isotropic model, where the singularities in the classical theory 
are point-like, the approach to singularity may be barrel, cigar or pancake type in anisotropic 
models. This could also potentially affect the pre-bounce attractor behavior. Further, in the case 
of anisotropic models, the phase of super-inflation becomes more non-trivial due to the bounces 
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in three different scale factors. It will be thus important to understand in future work the role of 
anisotropies on e-foldings obtained in the bounce regime of the power-law and assisted inflation 
scenarios in Bianchi models. The next step in understanding above results in a broader setting 
would require going beyond the homogeneity assumption. At the current stage of research, loop 
quantum effects as discussed in this manuscript are yet to be derived from loop quantum gravity 
and the way inhomogeneities affect the background dynamics is also not understood. However, 



important insights on connection of loop quantum cosmology and loop quantum gravity 60, l6ll |. 



and inclusion of inhomogeneities using spinfoam techniques have been gained 62], [63|. Though it 
leaves the discussion of affects of inhomogeneities on above results open, if the singularity resolution 
in loop quantum gravity occurs via a bounce, we expect some of the main results to remain 
unaffected, at least qualitatively. This expectation is based on the observation that a bounce 
of the homogeneous part of the metric would lead to a phase of super-inflation, which under the 
assumption that inhomogeneities remain small and dynamics remaining stable, would cause similar 
effects on power-law and assisted inflation scenarios. Whether or not this expectation holds, will 
be answered by future research on inclusion of inhomogeneities in this setting. 

Results obtained in this work enable us to rigorously understand the dynamics for exponential 
potentials in the presence of single and multiple fields, and bring forward some effects of the super- 
inflationary phase which were not noticed earlier with other forms of matter and potentials. In 
future work, it will be important to understand the effect that this phase and the increase in 
the number of e-foldings have on observational constraints on power-law and assisted inflation 
scenarios. The power-law inflation was recently briefly considered in the context of LQC [59(], 



however one would need to go beyond the underlying approximations in Ref. |59l ] to study in detail 
the implications of our results on the cosmological perturbations. It will be also interesting to 
extend this analysis by including different fluids, and study the effects of the scalar-fluid interaction 
on the dependence of the number of e-foldings during super-inflation on the steepness parameter 
and the equation of state. Finally, it will be instructive to analytically understand the approach 
of dynamical trajectories in LQC, particularly in the Planck regime, to the inflationary attractor 
in the future evolution. 
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